Magic angles and cross-hatching instabihty in hydrogel fracture 
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The full 2D analysis of roughness profiles of fracture surfaces resulting from quasi-static crack 
propagation in gelatin gels reveals an original behavior characterized by (i) strong anisotropy with 
maximum roughness at 1^-independent symmetry-preserving angles, (ii) a sub-critical instability 
leading, below a critical velocity, to a cross-hatched regime due to straight macrosteps drifting at 
the same magic angles and nucleated on crack-pinning network inhomogeneities. Step height values 
are determined by the width of the strain-hardened zone, governed by the elastic crack blunting 
characteristic of soft solids with breaking stresses much larger that low strain moduli. 

PACS numbers: 62.20.Mk, 83.80.Kn 



Over the past two decades, considerable effort has been 
devoted to characterizing and understanding the statisti- 
cal properties of the roughness of fracture surfaces in hn- 
ear elastic disordered materials. Investigation of a wide 
variety of systems, ranging from brittle sihca glass to 
ductile metallic alloys, has revealed the ubiquity of wide 
roughness spectra exhibiting self-affine characteristics on 
sizeable wavelength ranges. On the basis of analysis of 
height-height correlation functions in a single direction, 
E. Bouchaud et al. [1] first conjectured a fully universal 
behavior summed up into a single Hurst exponent. The 
discrepancies between the predictions of various subse- 
quent theoretical models have pointed toward the im- 
portance of also investigating possible anisotropics of the 
scaling properties • Work along this line indeed reveals 
different scaling exponents along the crack propagation 
direction and the (orthogonal) crack front one. From 
this, Ponson et al. [3| conclude that the self-affinity of the 
roughness is described by a Family- Vicsek scaling, hence 
by a set of two exponents. Within their new analysis, full 
universality is no longer the case since they evidence the 
existence of at least two classes of materials. 

However, on the basis of their theoretical work, Bouch- 
binder et al. Q have recently raised several new issues. 
In particular, they point to the necessity of a full 2D 
analysis of the correlations and find, when reexamining 
some experimental data, that a host of up to now uniden- 
tified exponents appear. They insist on the need to focus 
on the effects of departures from linear elasticity in the 
fracture process. 

In this spirit we report here on the nature of sur- 
face morphologies resulting from the fracture of gelatin, 
a highly compliant thermally reversible hydrogel, in the 
strongly subsonic regime. We show that these surfaces 
present very striking anisotropic features. Namely, the 
rms roughness R{0) , measured along a direction at angle 
from the propagation one Ox, exhibits two symmetry- 
preserving maxima for 9 ~ iOm- For gels with a fixed 
gelatin concentration, this "magic angle" 9m is constant 
over more than one decade of both crack velocity V 
and solvent viscosity rj. Moreover, below a critical, rj- 



dependent, velocity Vc a new regime develops, charac- 
terized by a cross-hatched (CH) morphology analogous 
to those observed previously on a covalently cross-linked 
hydrogel 0, @| and on various elastomers, swollen or not 
[7|. It is due to straight "macroscopic" steps, here of 
heights a few 100 fim, emerging from the above-described 
coexistent micrometric roughness. These steps are ori- 
ented at the same magic angles ±6*™ as those revealed 
by the microroughness. Their growth results from a dy- 
namical instability which we show to be subcritical: for 
V > Vc, while the inhomogeneities of the gel network are 
not strong enough to trigger step nucleation, we are able 
to force it with the help of externally applied local stress 
perturbations. 

These results lead us to propose the existence of a new 
class of morphologies, characteristic of highly compliant 
random polymer networks with breaking stresses larger 
than their small strain Young modulus, in which fracture 
implies elastic crack blunting [qI and strong elastic non- 
linearities at work in the near-tip region. 

Experimental — Our fracture surfaces result from 
mode-I cracks propagating along the mid-plane of long, 
thin gelatin slabs (Fig. l.a.), as described in detail in 
reference • A displacement A of the rigid grips is im- 
posed along Oy to the notched plate. After an initial 
transient, the crack reaches a steady regime at a veloc- 
ity V which increases with the "energy release rate" Q 
imposed by the opening A. We found that, in this slow, 
subsonic regime {V < 30 mm/s): 



e = 00 + rr]V 



(1) 



with Qq ~ 2.5 J.m~^, F ~ 10^, and were able to as- 
sign this strong ^-dependence to the fact that, in such 
a reversible gel, fracture occurs via scissionless chain 
pull-out, the high dissipation being due to viscous drag 

BE!. 

The results reported here were obtained on gels with 
a 5 wt.% content of gelatin in the solvent (pure water 
except when otherwise specified). The gel low-strain 
shear modulus fi = 3.5 kPa (transverse sound speed 
^ 1.9 m/s), is due to the entropic elasticity of the ran- 
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FIG. 1: Roughness characteristics for a 5 wt% gelatin/water gel for crack velocity V — 13.5 mm/s. (a) Experimental geometry. 
Sample dimensions -.11 — 30 mm, e = 10 mm, length L — 300 mm. (b) Height-height correlation functions Ahx{Ax) (lower 
curve) and Ah^iAz) (upper curve). Inset : a typical profile along Oz. (c) 3 x 3 mm^ map of heights h{x, z). The gray scale 
ranges over 27 /im. The oblique hnes correspond to peak positions in fig. l.d. (d) Anisotropic squared roughness E? vs. angle 
6 measured from crack propagation direction. The dark line is the best fit with a sum of two gaussians with an off'set baseline. 



tip velocity (mm.s^^) 


1.16 


5.9 


13.5 


23.3 


crack opening A (mm) 


5 


7 


9 


11 


rms roughness (/im) 


2.5 


3.1 


4.1 


5.3 


angle 9^ (°) 


18 


17 


17 


22 



TABLE I: 



dom polymer network, of average mesh size ^ ~ 10 nm. 
Note that, owing to this huge comphance, even very slow 
cracks {V ~ 100 /im/s) demand large openings A (Table 
I), typically > 3 mm. Post-mortem replicas of crack sur- 
faces, obtained by UV curing of a thin layer of glue, are 
scanned with the help of a mechanical profilometer (sty- 
lus tip radius 2 ^va, out of plane resolution 0.1 nm). Our 
maps h{x^ z) of heights above mean plane correspond to 
3x3 mm^ scans. Each one consists of a series of 600 
lines parallel to the front direction Oz, each line contain- 
ing 1500 data points. 

The microrough morphology — For V > 300 /im/s, 
the fracture surfaces are flat on the macroscale, yet 
clearly not mirror-like. As expected with such a ma- 
terial, elastic modulus as well as toughness fluctuations 
due to the randomness of the network result in in- and 
out-of-plane excursions of the crack front line, the profile 
roughness being the imprint of the latter ones. A typ- 
ical map is shown in figure l.b. For this crack velocity 
V — 13.5 mrn/s, we find that the global rms roughness 
R = 4.1 /xm [13|. Following Ponson et al. 3j, we first 
characterize height correlations along the front (Oz) i.e. 

1 /2 

we compute : Ahz{Az) — {[h{z + Az , x) — h{z , x)]) J ^ 
and its counterpart Ah^iAx) along the propagation di- 
rection (Ox). As indicated by figure l.b. the microrough- 
ness spectrum is anisotropic and of the broad band type. 
However, at variance with the findings of a simple 
power law fit (see dashed lines) is, here, clearly uncon- 
vincing. The discrepancy leads us, following Bouchbinder 



et al. [4^ , to try and exploit more extensively the full 2D 
information. Indeed, naked-eye inspection under grazing 
incidence, as well as the aspect of height maps (Fig. I.e.), 
strongly suggest the presence of preferred oblique direc- 
tions. This suspicion we have confirmed by computing 
the angle-dependent squared roughness B?{9), obtained 
by integrating the power spectrum |/i(q)p over an angu- 
lar sector dO = 4° about angle 9 + Tr/2 (in q-space) and 
in the wavelength window 2 /im-1 mm. We find in all 
cases that B?{0) exhibits two marked symmetric peaks 
at 6* = i^ni (see figure l.d. for which Om = 22°). More- 
over (Table I), Om does not show any significant variation 
over the whole V-range iVmaxlVmin 20). In addition, 
a 10-fold increase of solvent viscosity 77, obtained by using 
a 40/60 water/glycerol mixture, does not either affect its 
value. 

Finally (see Table I), the full 2D- averaged rms rough- 
ness R clearly increases with the imposed crack opening 
A, hence with V . 

The crosshatched ( CH) morphology — Starting from 
the microrough regime, we perform a step-by-step de- 
crease of A, thus slowing down crack propagation. Be- 
low a critical velocity Vc = 350 ± 20/im/s we observe on 
the fracture surfaces the emergence, at apparently ran- 
dom locations across the slab width, of oblique straight 
line defects (Fig. 2. a.). Upon further V decrease, these 
lines proliferate into a scale-like pattern, akin to the CH 
morphologies found by Tanaka et al. [5] on a covalently 
cross-linked gel and by Gent et al. 7] on elastomers. 

Note that, since these observations pertain to the qua- 
sistatic regime, front wave propagation [8,] is irrelevant 
here. 

Profile analysis (Fig. 2.b.) shows that these lines corre- 
spond to steps which emerge via self-amplification of the 
microroughness, with which they coexist. Their height 
saturates, after a traveling distance of several millime- 
ters, to values on the order of 200 /im. No change of 
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crystals, mortar, wood and sandstone Namely, none 
of these materials has been reported to exhibit (i) prefer- 
ential anisotropy directions (ii) a CH regime (iii) a notice- 
able increase of i? with V in the quasistatic regime. On 
the other hand, these a priori puzzling features are not 
restricted to gelatin gels. Indeed, as already mentioned, 
the CH regime has been documented for a chemical hy- 
drogel and for elastomers. Moreover, using figure 5 of 
reference Q, one easily checks that, in their chemical gel 
as well, the CH drift angle is constant {dm — 43°) over 
a 40-fold increase of y up to a critical Vc where the CH 
regime disappears. 




FIG. 2: Emergence of the CH regime, (a) A typical CH 
pattern close below Vc- The crack propagates upward, (b) 3D 
profile (400x400 fim^, 15 fim peak to peak) the emergence of 
a macrostep (arrow on Fig. 2. a.) out of the microroughness. 
(c) Critical velocity 14 vs. inverse viscosity for 5% gelatin gels 
in water/glycerol solvents. The line is the best linear fit. 



their drift angle is discernable as they grow, so that they 
propagate at the magic angle 9m of the microroughness. 
We find it to remain ^-independent over the full V < Vc 
explored range. 

The large A values make it possible to image the crack 
surfaces in the tip region. Figures 3.a-d depict the emer- 
gence and growth of a defect, in a frame attached to the 
moving crack front. It is seen to drift along the front 
direction Oz at a constant velocity Vd, as expected from 
the straightness of CH lines. So, as we have checked 
from e.g. Fig. 3.b-d., Vd — Vta.n9m- Fig. 3.e cor- 
responds to the fully developed defect, schematized on 
Fig. 3.f. It shows that the crack front has become dis- 
continuous, and composed of two semi-cracks AqA, BqB. 
As it moves along Ox, AqA creates two "lips" a/, a//. 
In the vicinity of its tip, the second semi-crack is de- 
layed with respect to AqA, so that its end point B lies 
on the an lip. The defect topology thus involves a ma- 
terial overhang ABCjDj. Note that this topology is the 
equivalent, for a widely open crack, of that identified in 
and 01 . As the front proceeds, the fold AiC'iDjBj 
ultimately collapses, leaving two complementary steps on 
the crack faces. 

Finally, in order to shed light on the role of solvent vis- 
cosity 77, we have studied gels made with water-glycerol 
mixtures. We find that the CH morphology persists but 
that the critical velocity Vc scales as I/77 (Fig. 2.c). How- 
ever, the magic angle 9m is 77-independent. 

Discussion — The above phenomenology presents 
several striking discrepancies with the known morpholo- 
gies, pertaining to metallic alloys, silica glasses, quasi- 




FIG. 3: Emergence and growth of a CH defect as viewed at 
60° from Ox in the Oxz plane. White bar length : 200 fim 
(a)-(d) Snapshots at constant time intervals At = 0.4s in the 
moving tip frame. The defect (bright line) grows and drifts 
at constant velocity, (e) Fully developed defect. The wite 
dashed lines mark the discontinuous front line, (f) Sketch of 
the defect topology. 

The question then arises of which material features 
may be responsible for this new class of behavior. In view 
of the dependence of Vc on solvent viscosity, one might 
be tempted to invoke poroelasticity. However Gent et al. 
have observed CH morphologies as well with dry elas- 
tomers as with a network swollen (ration 2.6) with oil. 
Hence the presence of solvent is not necessary for this 
regime to exist. 

This apparent paradox has led us to investigate the 
nature of the transition at Vc- The fact that (i) the mi- 
croroughness does not show any trend towards divergence 
when approaching Vc from above, (ii) the height 5 of CH 
steps does not vanish as ^ ^ V~ , prompts that the CH 
instability is discontinuous, i.e. involves nucleation. If so, 
it should be possible to trigger CH steps above Vc with 
the help of large enough external perturbations. Indeed, 
in a large range of V > Vc, inducing a mixed (II-I-III) 
mode stress perturbation by pinching the gel slab off mid- 
plane, slightly ahead of the tip, results in the nucleation 
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of a step that propagates through the sample. 

So, Vc appears as a threshold below which the struc- 
tural noise due the randomness of the gel network is not 
strong enough for the step nucleation barrier to be over- 
come. This, together with the fact that rjVc is a constant 
for a given network structure, suggests the following ten- 
tative picture: when the crack front meets a tougher zone 
(due to a local decrease of meshsize £_) of spatial ex- 
tent > dart, with dact ~ 100 nm the size of the process 
zone[10|, it slows down locally, generating a wall which 
stretches as the crack lips open. We postulate that step 
nucleation demands that the tough zone pins the front 
down to a full stop, i.e. that the front velocity fluctua- 
tion SV = —Vr- According to the fracture model devel- 
oped in refs. [13], [HI, both Qq and F in eq.dl]) scale as 
In the experiment, Q is kept constant via the con- 
stancy of A. So, -qSV = (^/r)(5fVC^- From this, we can 
estimate the maximum relative amplitude of mesh 
size fluctuations on spatial scales > dact m our system 
as Cm « r/20or]Vc « 7%. The smallness of em justifies 
our approximating Q by Qo in this estimate and explains 
the quasi-constancy of rjVc- Such an order of magnitude 
is statically plausible, and at the same time compatible 
with the sizeable large-angle light scattering exhibited by 
our gels. 

On the basis of this picture we interpret step nucleation 
as resulting from the pinning of a small part of the crack 
line. Due to local asymmetries, one of the resulting semi- 
cracks takes the lead, the other one getting delayed while 
skirting around the tough zone via an excursion towards 
the easier lip, thus nucleating the overhang. As cracking 
proceeds its endpoint B is gradually advected away by 
the opening of the AqA crack (Fig. 3.f), until it gets far 
enough from AqA for tensile stresses along an to become 
negligible. 

Now, our gelatin gels share with all the other soft 
solids exhibiting CH two mechanical characteristics: (i) 
the breaking stress is much larger than the low strain 
Young modulus E (ii) they strain-harden considerably 
at large deformation levels where polymer segments are 
stretched quasi-taut [13]. Hui et al. [9| have shown that 
(i) entails elastic crack blunting while, due to (ii) stress 
concentration is restricted to a strain-hardened region of 
thickness ^ Q/E along Oy. This allows us to predict that 
outward advection of B stops when the distance from B 
to AqA, i.e. the step height i5 becomes of order Qo/E. 
With Go K 2.5 J.m-2, £; « 10 kPa, we get 5 « 150 /im, 
in excellent agreement with the experimental value. For 
slow cracks in elastomers as studied in [?!, C/o and the 
relaxed E are respectively of order a few 10 J.m~^ and 
a few MPa, leading to step heights 5 in the 10 /im range. 
This agrees with Gent's findings, lending further support 
to our picture. 

Finally, let us reconsider the defect shape as depicted 
on Fig. S.f. In its end region close to A (i.e. for za < 
z < zb), the lower semi-crack connects the stretched lip 



aj with the ceiling of the fold, in which stresses are, to 
a large extent, relaxed by the BqB cut. This is consis- 
tent with the fact that, as cracking proceeds, BqB which 
cuts into the unscreened an lip, elongates {B drifts) to- 
wards z > while AqA regresses. This screening effect 
of the fold geometry explains that the defect drifts along 
the front in a direction defined by the initial symmetry- 
breaking choice which has led B^ towards one or the other 
crack lip. 

In the context of the above description, the existence 
of a preferred orientation in the micro-roughness could 
be attributed to aborted attempts at fold development. 
However, why 6m is independent of both V and the am- 
plitude of out-of-plane excursions remains at this stage a 
fully open question. 

In conclusion we propose that the phenomenology 
brought to light and analyzed here is characteristic, 
beyond the case of gelatin gels, of a new class of 
morphologies of fracture surfaces. The scenario which 
we propose to describe the cross-hatching instability 
implies that this class is characteristic of soft, tough, 
elastic materials which exhibit elastic crack-blunting and 
strong, strain-hardening, elastic non-linearities. 
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